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SUMMARY 


Development and evaluation of the Rockwell NASTRAN four-node quadrilateral 
(QUAD4) elanent is presented. The element derivation utilizes bilinear 
isoparametric techniques both for membrane and bending characteristics. The 
QUAD4 element coordinate system, membrane properties, lumped mass matrix, and 
treatment of warping are based upon the COSMIC/ NASTRAN QDMEM1 element while the 
bending characteristics are based upon a paper by T. J. R. Hughes. The effects 
of warping on the bending stiffness, consistent mass, and geometric stiffness 
are based upon a paper by R. H. MacNeal . Numerical integration is accomplished 
by Gaussian quadrature on a 2 x 2 grid. Practical user support features include 
valuable element thickness, thermal analysis and layered composite material 
definitions . 


INTRODUCTION 


Rockwell NASTRAN is the NASA/COSMIC released NASTRAN with Rockwell 
developed technical and efficiency enhancements incorporated. A total of nine 
Rockwell divisions fund the NASTRAN Group Service activities which include user 
consultation, development, maintenance, and validation of the production 
program. Rockwell NASTRAN is installed on IBM and CDC computing systems at 
three geograph. cal locations. The program is being used by the participating 
divisions which are located in California, Oklahoma, Ohio, Michigan and 
Pennsylvania. 
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The Rockwell QUAC4 has been developed in order to provide our users with a 
state-of-the-art general quadrilateral eleme.it. The improved efficiency and 
greater accuracy provided by this element eliminate the need of any of the other 
COSMIC/NASTRAN quadrilateral elements. Practical user support features 
incorporated in the development include varying element thickness, thermal 
strains, and laminated composite material inputs. The element derivation 
utilizes bilinear isoparametric techniques both for membrane and bending 
characteristics with numerical integration beirg accomplished by Gaussian 
quadrature on a 2 x 2 grid. 

The QUAD4 element coordinate system, membrane properties, lumped mass 
matrix and treatment of warping are based upon the COSMIC/NASTRAN QDMEM1 element 
while the bending properties are based upon a recent paper by T. J. R. Hughes 
(ref. 1). The effects of warping on the bending stiffness, consistent mass, and 
geometric stiffness are based upon a paper by R. H. MacNeal (ref. 2). The 
theory adopted from reference 1 appears to minimize or preclude some of the 
complications alluded to in reference 2. In particular, no special local 
Cartesian system or selective integration procedure is required to achieve a 
reasonably good element behavior. 

General theoretical background of the element stiffness matrix is presented 
in equations 1 through 35 of the theoretical background section. Derivation of 
the equivalent thermal aDplied load vector is presented in equations 38 through 
41. 


The evaluation of element test results as proposed by reference 3 are 
presented ir. table 1. The test results for static analysis of various 
structures, mechanical loadings, and thermal analysis are presented in tables 2 
through 8. The results for the real eigenvalue test case are presented in table 
9. The transverse central deflection computed for the three composite material 
test cases using the MSC/QUAD4 and the Rockwe 1 1 /QUAD4 element is presented in 
table 10. 


THEORETICAL BACKGROUND 


The relationship between forces and strains (including thermal terms) is 

aesc-ibed by the following matrix where the vectors (eS and {x*} are thermal 
generated strains and curvatures, respectively. 
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The limits on the integration are from the bottom surface to the top 
surface of the plate. The matrix of material moduli, [Ge] , has the following 
form for orthotropic materials: 
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Here, vj ~ v 2 E 1 ’ is re( l uired the matrix of elastic moduli be 
symmetric. The [G^] is a 2 x 2 matrix of elastic coefficients for transverse 
shear. H s , the effective thickness for transverse shear, has a default value of 
H $ / H = 5/6, which is the correct value for a homogeneous plate with an 
actual membrane thickness of H. 

Figure 1 depicts a plate composed of the eight laminas. For this case. A, B and 
D are defined as follows: 


A = 

' £ G e ! h k ~ h k-i) 

(13) 
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Let A and N denote the area and shape functions, respectively, of an element, 

d 

where n is the number of element nodes. For the case of a homogeneous, 
isotropic, linearly elastic plate of thickness H, the element stiffness matrix, 
K e , may be defined as follows. 
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The formulation of the element stiffness matrix follows the procedure 
defined in reference 4 and 5. Then R b, s can be written in the following form: 


R 


b 

a 


0 0 N a ': 

0 N al 0 

0 N a >2 N a ’: 


1 < a < n (21) 


The shear stiffness is obtained by the technique mentioned in reference 1. 
The detailed procedures are discussed next. 

Geometric and kinematic data are defined in figure 2, and the direction 
vectors have unit length (e.g. !l e njj = Let w a ^ © a denote the 

transverse displacement and rotation vector, respectively, associated with node 
a. Throughout, a subscript b will equal a+1 modulo 4. 

The definition of the element shear strains may be described in the 
followir steps. 

(1) For each element side, define a shear strain component at the 
midpoint, in a direction parallel to tne side. 

g a = H-W a )/ eal-'^V^a 1 *Z (22) 
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(2) For each node, define a shear strain vector. (See figure 3 geometric 
interpretation of this process.) 


V b = 7 bi e bl +Y b2 e b: 

(23) 

Ybi =,10, b' |1( gbi'gb2“ b) 

(24) 

yb 2 =,1 ’*b )- ^^'gbiO'b 1 

(25) 

«b =e bf e b: 

(26) 

gbi = gb 

(27) 

®b2 a 

(28) 


(3) Interpolate the nodal values by way of the bilinear shape functions 

'V s ) 

r-i.Va (29) 

3— i 

For the transverse shear strain interpolations derived in the previous section, 
R s takes on the following form: 


R b 

=Rbi R b 2 R b 3 15b - 4 

(30) 
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The element stress resultants may be obtained from the following relations: 
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where 


element displacement vector 


Finally, thermal expansion is represented by a vector of thermal strains 
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(T - Tn) = 
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( 38 ) 


where a ~ thermal expansion coefficients 

T = Temperature at any point in the element 
T q = reference temperature of the material 

An equivalent elastic state of stress that will produce the same thermal 
stress is 


M- icit' 1 ! 


( 39 ) 
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An equivalent set of generalized loads P applied to grid points of the 
element is obtained by 


P 



(40) 


The equivalent thermal moment vector is defined as 


M t 



T' 


z dz 


where T' is the thermal gradient at a cross-section of the plate. 


(41) 


numerical examples 


The test problems have been selected from reference 3. The elements tested 
included the C0SMIC/QUAD2, the MSC/QUAD4 and the Rockwell/QUAD4. The test runs 
for QUAD2 and QUAD4 were performed on an IBM 3081 computer at the Rockwell 
Western Computing Center while the MSC/Q*'A04 result were obtained by utilizing 
version 63 of MSC/NASTRAN on the Rockwell Scientific Computing Center CDC/CYBER 
equipment. 

The grading system for finite elements proposed by reference 3 is: 

Grade Range 

A 2% > Error 

B 10% > Error > 2% 

C 20% i Error > 10% 

D 50% 2 Error > 20% 

F Error > 50% 

The structures analyzed to evaluate the test elements included a patch test 
plate (figure 4), a straight cantilever beam (figure 5), a curved cantilever 
beam (figure 6), a rectangular plate with different aspect ratios (figure 7), a 
Scordelis-Lo roof (figure 8), and a simply supported plate (figure 9) for normal 
modes and layered composite analysis. 
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Table 1 Drpcertts Ine summary of grading results for the tested elements. 
Thp result for each of the individual test cases are reported in tables 2 
through 10. The patch test results presented in Table 2 are reported in the 
form of percentage error of the computed stresses. The results reported in 
tables 3 through 7 are shown in normalized form where the computed displacement 
data has been divided by the theoretical value. The most disturbing failure of 
the QUAD2 element is its inability to get a passing grade for the straight beam 
in-plane shear and twist cases. QUAD2 also failed in the curved beam and 
Scorde'lis-Lo roof problems. Neither of the QUA04's or the QUAD2 could achieve a 
passing grade for the straight beam in-plane shear with trapezoidal shaped 
elements. In general, our published results agree, but there are some 
differences from those reported in reference 3. In particular, the results of 
the twist case for all element configurations of the straight cantilevered beam 
problem do not agree with the results presented in reference 3. We believe that 
this was due to a problem with version 63 of MSC/NASTRAN as installed on our CDC 
equipment at the time we were making our test case runs. 


CONCLUSION 

In this paper, we have examined the behavior of the new four-node 
quadrilateral element implemented in Rockwell NASTRAN. The element has been 
shown to behave well for a variety of pla.e problems and has retained simplicity 
in the formulation. The formulation enabled straightforward generation of a 
linear triangular bending element, which has also been successfully implemented 
in Rockwell /NASTRAN. 
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Table 1 Summary of Test Results 


Test 

Table 

Element 

Shape 

"K 17“ 
QUAD4 

QUAD2 

user 

QUAD4 

1 Patch Test, Membrane 

2(a) 

Irregular 

A 

A 

A 

2 Patch Test, Bending 

2(b) 

Irregular 

A 

A 

A 

3 Straight Beam, Extension 

3(a,b>c) All 

A 

A 

A 

4 Straight Beam, In-Plane Shear 

3(a) 

Regular 

B 

F 

B 

5 Straight Beam, In- Plane Shear 

3(0 

Irregular 

D 

F 

D 

6 Straight Beam, Out-of-Plane Shear 

3(b,c) 

Regular 

A 

B 

A 

7 Straight Beam, Out-of-Plane Shear 

3(b,c) 

Irregular 

B 

B 

B 

8 Straight Beam, TVist 

3(a,b,c) All 

B 

D 

D 

9 Curved Beam, In-Plane Shear 

4 

Regular 

C 

F 

C 

10 Curved Beam, Out-of-Plane Shear 

4 

Regular 

B 

D 

C 

11 Rectangular Plate (N=4) 

5.6(a) 

Regular 

A 

A 

A 

12 Scordelis-Lo Roof (N=4) 

7 

Regular 

A 

D 

A 


Failed Test Grade (D's and F's) 

Table 2 Patch Test Results (Figure 4~J 
Max. Errors (%) of Stress 


(a) Membrane Plate 



RI/QUAD4 

quAD2 

MSC/QUAD4 

(X\ - Oy 

4.2 

4.2 

0.0 

T X V 

/ 

1.0 

1.0 

0.0 

(b) Braiding Plate 




in x s Hly 

4.2 

4.2 

0.0 

& 

X 

0.9 

0.9 

0.0 

y = (Ty 

4.2 

4.2 

0.0 

T *y 

1.0 

1.0 

0.3 


57 



Table 3. Results for Straight Cantilever Beam (Fig. 5) 
Normalized Tip Displacement in Direction of Load 


Tip Loading 
Direction 

RI/QUAD4 

QUAD?. 

MSC/QUAD4 

(a) 

Rectangular Elements 



Extension 

0.996 

0.992 

0.996 

In- Plane Shear 

0.904 

0.032 

0.904 

Out-of-Plane Shear 

0.980 

0.971 

0.986 

TVist 

0.941 

0.567 

0.702 

(b) 

Trapezoidal Elements 



Extension 

0.996 

0.993 

0.996 

In- Plane Shear 

0.071 

0.016 

0.071 

Out -of -Plane Shear 

0.964 

0.963 

0.958 

Twist 

0.884 

0.605 

0.705 

(c) 

Parallelogram Elements 



Extension 

0.996 

0.992 

0.996 

In- Plane Shear 

0.808 

0.144 

0.795 

Out-of-Plane Shear 

0.978 

0.961 

0.977 

TWist 

0.849 

0.615 

0.705 

Table 4 Results for Curved Beam (Fig. 6) 
Normalized Tip Displacement in Direction of Load 


Tin Loading Direction RI/QUAD4 

qUAD2 

MSC/QUAD4 

In- Plane Vertical 

0.835 

0.025 

0.835 

Out-of-Plane 

0.956 

0.597 

0.868 

Table 5 Results for Rectangular Plate Simple Supports (Fig. 7) 


with Concentrated Load 


Normalized Transverse Deflection at Center 



(a) Aspect Ratio = 1.0 



Mesh Size(N)* 

RI/QUAD4 

QUAD2 

MSC/QUAD4 

2 

0.992 

1.035 

0.960 

4 

0.995 

1.011 

1.017 

8 

1.033 

1.083 

1.045 


(b) Aspect Ratio =5.0 



2 

0.844 

0.493 

0.870 

4 

0.928 

0.685 

0.962 

8 

0.986 

0.845 

1.005 


* only one quadrant is discretized 
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Table 6 Results for Rectangular Plate Clamped Supports 
(Figure 7) With a Uniform Load 


Normalized Lateral Deflection at Center 


(a) Aspect Ratio = 1.0 


Mesh Size (N) 

RI/QUAD4 

QUAD2 

MSC/QUAD4 

2 

0.961 

1.024 

1.008 

4 

0.993 

1.019 

1.032 

8 

1.016 

1.057 

1.040 


(b) Aspect Ratio = 5.0 



2 

1.124 

0.873 

1.314 

4 

0.962 

1.001 

1.016 

8 

1.002 

1.019 

1.0’" 


Table 7 Results For Scordelis-Lo Roof (Figure 8) 


Normalized Vertical Deflection at Midpor^ of Free Edge 


Mesh Size (N) 

RI/QUAD4 

qUAD2 

MSC/QUAD4 

2 

1.309 

0.881 

1.313 

4 

1.017 

0.690 

1.021 



Table v Conpar ison of Analytical, QUAD4, 
NASTRAN D04O 1-11-1 (Reference 6) 

Max. 

and QUAD1 

Category 

Analytical 

CQUAD1 

CQUAD4 


-1 

-1 

-1 

Displacement 

6.2898x10 

6.2895x10 

6.317195x10 

2 

2 

2 

Moment m v 

1.4770x10 

1.4888x10 

1.4832200x10 

/ 

3 

3 

3 

Stress Ty 

7.764618x10 

7.792977x10 

7.779586x10 
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Table 9 Natural Frequency Comparison, cps or 
N^STRAN/DEMO 3-1-2 (4x4) 


Mode No. 

Theoretical 

RI/QUAD4 

qUADl 

1 

0.9069 

0.8823 

0.9021 

2 

2.2672 

2.3376 

2.2837 

3 

4. 5345 

4.3515 

4.7179 


Table 10 Transverse Central Deflection of Simply Supported 
Conposite Square Plate Under a Uniform Pressure (Figure 9) 


of Plies 

Type of Lcminate 

RI/QUAD4(IBM) 


2 

90 °/ 0 ° 

-3 

5.63410x10 

-3 

.,.589612x10 

3 

0°/90°/0° 

-3 

5.55896x10 

-3 

5.423785x10 

4 

90°/0°/90°/0 c 

-3 

5.58961x10 

-3 

5.666982x10 
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a = .12; b = .24; t = .001 
fc = 1.0 x 10 6 ; •- = 0.25 

Location of Inner Nodes: 

x y 

1 .04 .02 

2 .18 .05 

3 .16 .08 

4 .08 .08 

Boundary Conditions: 

(a) Membrane 

u = 10° (x + y/2) 
v = 10° (y + x/2) 

(b) Bending 

k = 10 (x“+xy+y')/2 

• x - io- (y+x/2) 

0 y = 10' 3 (-x-y/2) 


Figure 4. Patch Test for Plates 
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T 


a 


(a) Regular Shape Elements 



(b) Trapezoidal Shape Elements 



(c) Parallelogram Shape Elements 

Length = 6.0; Heig.it = 0.2; Thickness =0.1 
E = 1.0 x 10 7 ; *- = 0.3; Mesh = 6x1 
Loading: Unit forces at free end 



Figure 5. Straight Cantilever Beam 
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Inner radius = 4.12; Outer radius =4.52 
Arc = 90°; Thickness =0.1; ► = 0.25 
E = 1.0 x 10 7 ; Mesh = 6x1 
Loading: Unit forces at tip 





Figure 6. curved Beam 




a = 2.0; b = 2.0 or 10.0; ► = 0.5 , 

Thickness =0.001; E = 1.7472 x 10 ' 
Boundaries = simply supported or clamped 
Mesh = N x N (on 1/4 of plate) 4 
Loading: Uniform pressure o = 10 or 
Central load p = 4.0 x 10 4 


Figure 7. Rectangular Plate 



Radius = 25.0; Length = 50.0; Thickness =0.25 
* = 0.0; Loading =90.0 per unit area in -z direction 
E = '.32 x 10®; U = U = 0 on curved edges 

X 2 

Mesh = N x N on shaded area 


Figure 8. Scordelis-Lo Roof 
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ORIGINAL PAGE 
OF POOR QUALITY 



i = 10.0; T = 0.2; * = .25 

E 1 = 20. x 10 6 ; E, = .5 x 10 6 ; G = .25 x 10 6 

Loading Condition: 0.5 psi uniform pressure 

Case 1 : 2 piles, material angle of fiber 90°/0 c 

Case 2 : 3 piles, " " 0°/90 o /0° 

Case 5 : 4 piles, " " 90 o /0°/90 o /0° 


Figure 9. Simply Supported Square Plate for a Layered Structure 
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